Improved Limits on Spin-Dependent WIMP-Proton Interactions from a Two Liter 

CF 3 I Bubble Chamber 
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Data from the operation of a bubble chamber filled with 3.5 kg of CF3I in a shallow underground 
site are reported. An analysis of ultrasound signals accompanying bubble nucleations confirms that 
alpha decays generate a significantly louder acoustic emission than single nuclear recoils, leading to 
an efficient background discrimination. Three dark matter candidate events were observed during an 
effective exposure of 28.1 kg-day, consistent with a neutron background. This observation provides 
strong direct detection constraints on WIMP-proton spin-dependent scattering for WIMP masses 
> 20 GeV/c 2 . 
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There is abundant evidence that ^85% of the matter 
in the Universe is cold, dark, and non-baryonic pj. To- 
gether, dark matter and dark energy are pillars in mod- 
ern cosmology which form a simple framework to un- 
derstand the detailed observations of baryonic structure, 
the Hubble diagram, the cosmic microwave background, 
and the abundances of light elements. However, present 
knowledge of the nature of these essential ingredients of 
cosmology is limited. A new weak-scale symmetry, such 
as the preservation of R-parity in supersymmetry, would 
predict a massive particle with properties that conform 
with the current knowledge of dark matter Q. If these 
Weakly Interacting Massive Particles (WIMPs) are the 
dark matter, then they may scatter off nuclei with enough 
energy, and at a high enough rate, to be detectable in the 
laboratory Q. 

The Chicagoland Observatory for Underground Parti- 
cle Physics (COUPP) employs the bubble chamber tech- 
nique to search for WIMP-nucleon elastic scattering Q. 
The bubble chamber is a powerful device to explore nu- 
clear recoils in a rare event search. If the chamber pres- 
sure and temperature are chosen appropriately, electron 
recoils from the abundant gamma-ray and beta-decay 
backgrounds simply do not nucleate bubbles [f|. Neu- 
trons from spontaneous hssion and (alpha, n) in materi- 
als at the experimental site can be moderated with low- 
Z shielding materials, and cosmogenic neutrons can be 
reduced to negligible levels with an appropriate overbur- 
den. Due to the threshold nature of the bubble nucle- 
ation process, alpha-decays are more of a background 
concern for dark matter searches with superheated liq- 
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uids than they are for technologies with event- by- event 
energy measurements. Recently, however, the PICASSO 
collaboration reported event- by- event alpha-decay iden- 
tification based on the acoustic emission from bubbles 
6]. This letter confirms the observation of alpha dis- 
crimination and utilizes this technique to set limits on 
spin-dependent WIMP-proton scattering. 

This letter reports results from a 3.5 kg CF3I bub- 
ble chamber operated from August 19th to December 
18th, 2009, in the MINOS near detector tunnel Q at 
the Fermi National Accelerator Laboratory. The bubble 
chamber was located 3 meters off the axis of the NuMI 
neutrino beam Q, which operated with 10 /is pulses, 
typically every 2.5 seconds. The short duty factor al- 
lowed for a calibration source of efficiently tagged, beam- 
induced, fast neutrons. Additional calibration was pro- 
vided using a switchable americium-beryllium neutron 
source (sAmBe) @. 

The 300 ft overburden of the NuMI site was sufficient 
to shield out the hadronic component of the cosmic ray 
flux and to significantly attenuate the cosmic ray muon 
flux. Residual cosmic ray muons passing near the bubble 
chamber were tagged using a 1000 gallon Bicron 517L 
liquid scintillator counter equipped with 19 RCA-2425 
photomultiplier tubes which surrounded the experiment 
on the sides and above. The liquid scintillator, in con- 
junction with passive high-density polyethylene shield- 
ing below, provided a moderator for neutrons originating 
from natural radioactivity in the cavern rocks. 

The bubble chamber consisted of a 150 mm diameter 
3-liter synthetic fused silica [8] bell jar sealed to a flexible 
stainless steel bellows and immersed in propylene glycol 
within a stainless steel pressure vessel. The propylene 
glycol, which served as the hydraulic fluid to manage 
the inner pressure of the bubble chamber, was driven by 
an external pressure control unit. The flexible bellows 



served to ensure that the contents of the bell jar were 
at the same pressure as the hydraulic fluid. The bell 
jar contained 3.5 kg of CF3I topped with water. The 
water provided a buffer to ensure that the CF3I contacted 
only the smooth silica and not the rough stainless steel 
surfaces. 

The thermodynamic conditions of the chamber were 
monitored with two temperature sensors mounted on 
the bellows flanges and with three pressure transduc- 
ers. Four lead zirconate piezoelectric acoustic transduc- 
ers were epoxied to the exterior of the bell jar to record 
the acoustic emissions from bubble nucleations, the audi- 
ble "plink" used to trigger the flash lamps in early bubble 
chambers [9] . Two VGA resolution CMOS cameras were 
used to photograph the chamber with a 20-degree stereo 
angle at a rate of 100 frames per second. Frame-to-frame 
differences in the image data provided the primary trig- 
ger for the experiment, typically initiating a compression 
within 20 ms of a bubble nucleation. Stereo image data 
from the cameras were used to reconstruct the spatial 
coordinates of each bubble within the chamber. 

The chamber was operated at a temperature of 29.5°C 
and a pressure of 26.5 psia when sensitive to particle 
interactions, providing a Seitz model bubble nucleation 
threshold of 21 keV nuclear recoil energy [lOj. Although 
this threshold depends on a theoretical model, it has been 
benchmarked by comparing to the threshold for sensitiv- 
ity to alpha decays [5[. Reasonable uncertainties in the 
measurements of the thermodynamic state of the super- 
heated target, about 1 psi and 0.5°C, correspond to a 
15% shift in the threshold of the chamber and a 30% 
(15%) change in sensitivity for 50 (100) GeV WIMPs. 

Following the nucleation of a bubble, the chamber was 
compressed to 215 psia for 30 seconds (300 seconds ev- 
ery 10th event) to re-condense the CF3I vapor. During 
the compression period, data from the bubble event were 
logged, and the chamber was prepared for the next ex- 
pansion. To allow for equilibration of the bubble cham- 
ber fluid, the valid live time for the chamber began 30 
seconds after an expansion. From the outset, the per- 
formance of this chamber was a significant improvement 
over previous COUPP bubble chambers [5j. There was 
no evidence of radon beyond a few dozen events immedi- 
ately after the chamber fill. Vessel wall nucleations due 
to alpha decays, observed at a rate of 0.7 events cm~ 2 
day" 1 in natural quartz vessels, were at an unmeasurably 
low level in this synthetic fused silica vessel. 

The data presented in this letter confirm PICASSO's 
observation of a difference between the acoustic emissions 
from bubbles nucleated by nuclear recoils and bubbles nu- 
cleated by alpha decays [6j. The current understanding 
of the physics behind this discrimination is that the ul- 
trasonic acoustic emission from a growing bubble peaks 
when the bubble is larger than a nuclear recoil track, tens 
of nanometers, but smaller than an alpha track, tens of 
microns [ll(. At the time of peak acoustic emission, a nu- 
clear recoil track is a single bubble, while an alpha track 
may contain several microscopic bubbles making the al- 
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FIG. 1: The acoustic parameter, defined in the text, is a mea- 
sure of the loudness of the acoustic "plink" resulting from 
bubble nucleations. The red, dashed histogram in the main 
panel is from the nuclear recoil calibration datasets (sAmBe, 
NuMI, cosmic). The blue, solid histogram in the main panel 
is from all the untagged data, including data taken before 
the muon veto was commissioned. The inset histogram con- 
tains events from 222 Rn early in the run, which reproduce the 
second peak in the untagged dataset. The acceptance of the 
nuclear recoil cut, indicated with arrows, is 



pha event louder. Because the COUPP bubble chamber 
fluid is a homogeneous medium, it is expected that acous- 
tic alpha/nuclear recoil discrimination should be more 
distinct than in the aqueous gel emulsion of superheated 
droplets used by PICASSO. 

The acoustic transducer signals were digitized with a 
2.5 MHz sampling rate and recorded for 100 ms for each 
event. The sound of bubble nucleation showed a broad 
emission distinctly above background noise up to a fre- 
quency of 200 kHz. The analysis of these signals was 
based on an acoustic parameter (AP) which is a mea- 
sure of loudness. The AP is a frequency weighted acous- 
tic power density integral, corrected for sensor gain and 
bubble position 1 . The AP was scaled to have a value of 
unity at the peak observed in its distribution for bubbles 
induced by known neutron sources (the sAmBe source, 
the pulsed NuMI beam, or cosmic muons). The AP dis- 
tribution from these known neutron sources is shown in 



AP = J2j Gj J2 n C n (x) ^^" x / x psd?f. Cn(x) are corrections 
for the position dependence of the bubble acoustic emission, x 
is the position of the bubble calculated from the stereoscopic 
images. Gj is the gain factor for acoustic transducer j. / is 
frequency, n = 1 . . . 3 indicates three frequency bands (5-30 kHz, 
30-80 kHz, or 80-120 kHz), j identifies the acoustic sensor, / m j n 
and fmax are the boundaries of the frequency band, and psd 3 ^ 
is the power spectral density for the bin with center frequency / 
for sensor j. 
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Fig. Q] Also shown is the distribution of untagged data, 
mostly the ^300 kg-day exposure before the scintillator 
veto was installed, which showed a peak at AP = 1 simi- 
lar to the calibration data, with an additional component 
with 3 < AP < 5. The low AP peak in the untagged data 
is caused by nuclear recoils and the higher AP peak is 
understood to be caused by alpha decays. 

All data have been subject to a fiducial volume cut 
requiring the stereoscopically reconstructed bubble posi- 
tions to be at least 5 mm away from the wall. The ac- 
ceptance of this fiducial volume cut was 75% for nuclear 
recoils from the sAmBe and NuMI calibration samples. 
The definition of AP and of the AP cut were specified be- 
fore the veto was installed, and neither was tuned based 
on the WIMP search data sample. The nuclear recoil ac- 
ceptance of the AP cut was measured to be 88% in the 
fiducial volume using the NuMI coincident and sAmBe 
calibration events. 

There was no viable method of injecting alpha-emitting 
isotopes into this bubble chamber. On the contrary, ef- 
forts were focused on removing these from the active vol- 
ume, achieving a rate of ~ 0.7 alpha decays per kg day. 
A small amount of 222 Rn injected during the fill of the 
chamber was identifiable via its characteristic decay time. 
These events were selected by requiring a time difference 
between bubbles of < 3 minutes, corresponding to the 
sequence of alpha decays from 222 Rn decaying to 218 Po 
(Ti =3.1 min) which then emits a second alpha particle. 
This small sample of 10 known alpha decay events shown 
in the inset of Fig. [1] is consistent with the peak observed 
in the untagged data at 3 < AP < 5, confirming the 
expectation of a louder alpha acoustic emission. 

The camera images and the slow pressure rise asso- 
ciated with bubble expansion provided two methods to 
count macroscopic bubbles in an event. These bubble 
counting methods were used to produce Fig. [5J where 
the AP is plotted as a function of bubble number in- 
ferred from the pressure rise. The symbols indicate how 
many bubbles were inferred from camera images. WIMP 
interactions will nucleate only one bubble. Neutrons are 
the only significant background able to produce multi- 
ple macroscopic bubbles during an event. The AP scales 
with the number of macroscopic bubbles. Tails in the 
AP distribution for a given number of macroscopic bub- 
bles are all towards louder values. 2 This is consistent 
with the AP scaling linearly with acoustic energy, with 
a minimum acoustic energy produced by single bubbles. 

Previous COUPP calibration data taken at 30° C 
showed that neutron induced nuclear recoils with energy 
above the Seitz threshold had a 50% chance of nucleat- 
ing bubbles Q. A comparison of MCNP simulations of 



2 The high AP tails are a possible indication of inefficiency for 
the AP neutron selection, but not a background concern. The 
nature of the events in the tail of the acoustic parameter has not 
been studied. The rate is consistent with the expected rate of 
(n, alpha) reactions on flourine and carbon. 
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FIG. 2: The two methods of measuring the bubble multi- 
plicity are plotted for the full data set of 291 kg-day expo- 
sure. The first method, based on automated video reconstruc- 
tion software, searches for clusters of pixels changing between 
frames and is represented by the various symbols labeled in 
the legend. The second method of counting bubbles was to 
measure the dP/dt associated with an event. The abscissa 
is the pressure rise, scaled to the number of bubbles. The 
ordinate is the acoustic parameter AP (see text). 
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FIG. 3: Data from a 28.1 kg-day WIMP search, plotted as 
in Fig. [2] In the alpha region (large dashed box) 21 events 
with no evidence for muon activity were observed, with 3 
events having coincident muon activity (labeled vetoed). The 
smaller dashed box indicates the WIMP signal region, con- 
taining 3 candidate events and 47 events coincident with 
muons in the veto. The neutron double scatter region, in- 
dicated with a dashed circle, contains 8 vetoed events and 1 
unvetoed event, the latter indicating an irreducible neutron 
background contamination (see text). 
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FIG. 4: Improved limits on spin-dependent WIMP-proton 
elastic scattering from the data presented in this letter. A 
previous COUPP result [f| is shown for comparison. Direct 
detection limits from the PICASSO experiment [jig , cyan, 
and the KIMS experiment [l3l ]. green, are shown. Limits on 
neutralino annihilation in the Sun from the IceCube [TH, ma- 
genta, and Super Kamiokande [14j. black, neutrino observa- 
tories are also plotted. The indirect detection limits from 
the neutrino observations have additional dependence on the 
branching fractions of the annihilation products. The gold 
region indicates favored regions in cMSSM 16|. The blue, 
dashed-dotted line is the expected cross-section for "maver- 
ick" dark matter with Q.h 2 =0.1 [13 . 

single and multiple rates with the rates observed during 
sAmBe calibration is consistent with this 50% efficiency 
for the data presented here. The efficiency is expected to 
reach 100% at higher operating temperatures Q; unfor- 
tunately, a mechanical failure forced an early termination 
of the run, restricting data-taking to 30°C. 

The WIMP search data spanned November 19th to 
December 18th. They are shown in Fig. [3] along with 
the signal acceptance regions for various classes of event. 
There were 23 days of livctimc with the muon veto ac- 
tive. A number of data quality cuts were imposed to 
remove periods of poor detector performance. The most 
restrictive data quality cut was imposed to excise inter- 
mittent noise on the acoustic signals, removing 42% of 
the science run. An effective exposure of 28.1 kg-days 
remained after imposing all data quality cuts, the fidu- 
cial volume cut, and the AP nuclear recoil cut. In the 
signal acceptance region there were 3 WIMP candidate 
events and 47 events coincident with muon activity in 
the veto. In the alpha region there were 24 events, 3 of 
which were coincident with muon activity. These are sta- 
tistically consistent with the high-AP tail expected from 
the 47 cosmic ray events in the signal region. Of 9 events 



with two macroscopic bubbles, one was not coincident 
with the veto. This event is proof of an unvetoed neu- 
tron background component with sufficient rate to result 
in the observed 3 single-scatter WIMP candidate events. 

Assessments of event-by-event discrimination against 
alpha decays and of the sensitivity to WIMP-nucleon 
scattering are limited by the unvetoed neutron back- 
ground in this dataset. Interpreting the three events in 
the signal region as alpha decays results in a conserva- 
tive 90% C.L. upper limit on the binomial probability 
of an alpha decay registering in the nuclear recoil signal 
region of < 26%. At this operating pressure and tem- 
perature, an alpha particle will create bubble nucleation 
sites along its entire track, and there is clear evidence of a 
neutron background from the multiple scatter events, so 
these three events are likely not alpha decays. Therefore 
the presently derived alpha background rejection should 
be considered a conservative assessment for the potential 
of this technique. We expect an improved estimate from 
runs in a deeper underground site, where the residual 
neutron background should be absent. 

Interpreting the three events in the signal region as 
WIMP candidates results in a 90% Poisson upper limit 
of 6.7 for the mean of the signal. The resulting im- 
proved limits on spin-dependent WIMP-proton couplings 
are shown in Fig. HJ The spin-independent sensitivity 
that can be extracted from present data is comparable to 
that obtained by CDMS in another shallow underground 
facility [llj]. The calculations assume the standard halo 
parameterization fl9| . with pz> = 0.3 GeV c~ 2 cm~ 3 , 
v esc = 650 km/s, ve = 244 km/s, t>o = 230 km/s, and the 
spin-dependent couplings from the compilation in Tovey 
et a/.[20(. This result is consistent with a background 
from neutrons induced by residual cosmic radiation in 
the shallow site. 

In view of the ~ 10 -11 intrinsic rejection against min- 
imum ionizing backgrounds 5] and the acoustic alpha 
rejection demonstrated in this letter, a leading sensitiv- 
ity to both spin-dependent and -independent WIMP cou- 
plings can be expected from the operation of CF3I bub- 
ble chambers deep underground. A first exploration of 
the cMSSM spin-dependent parameter space [ll| of su- 
persymmetric dark matter candidates is expected from 
operation of this chamber in a deeper site. At the time 
of this writing, a 60 kg CF3I COUPP bubble chamber is 
being commissioned. 
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